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Chapter I 



PITRODUCTION 



1* Backgroiind of the Problem. 

From the moment a projectile leaves th- muzzle of 
the gun which launches it on its fli^t *»down range", it is 
subjected ,o a complicated variety of forces. The study 
oi the effects of these forces x^jon wh re the projectile 
goes, and upon how long it takes to get there, has for cen- 
turies fascinated mathematicians and physicists alike. 

Modern guns are- rifled, ai.d hence impart to the 
projectile a spinning motion which Ktt as its "stability" 
and thus influences its flight. But, even if one is 
willing to defer the tempting study of the gyroscopic 
properties of a spinning shell, on the ground oha^ th y 
are jecondary in importance to he mo ion o its center 
of mess, one finds tuat, to comp'.te a trajectory ot t: e 
raa V center, one must solve a di .ferential system involving 
a number of functions which depe d upon the place at which 
the proje-^ti'e happen to be as we 1 as upon th*j time when 
it happens to be the^e, not to menti n the spe* i nt which 
it happens to be moving. 

For ezamile, the air resistance whicr a bullet 
nccun*ers is dependent not only upo. the proj-ctiln’s 
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3pee< relative to tne eir bu clso upoa tl.e -.tio spheric 
d^’isity and temperet ‘"s, and th*.**e lettei quautli ies vary 
from place to pla* v and from Ins' ant to Instoat q ;5 the 
shell tx^versaa ^ta path* The alr-aueed of the projeotile 
is itaelf cieyf.nd:>iv upon the speed and direotlon of the 
wind» >nd v»ind8 are notoriously "oontrdr:^". 

At th^s wrltlngi not even with the aid of high- 
speed omputlng maohlnery has it bee. praotioable to obtain 
a oomplete tabulation of solutions of the differential equa- 
tions of motion* eovoring all possible sets of initial con- 
ditions, wind- struct urea, density and temperature distribu^ 
tions eto ?he evioe which ballistician^ have employed 
to overcome this difficulty is that of integrating the equa- 
tions under certain relatively simple •standard*, or •nor- 
ma!^* conditions (such as zero wind, a n n-rotating earth, 
eto,], and then conqnitlngi a 1 Inear approximation to the 
actual effect (on range, on tlae-of- 'light, or on some 
Other trajectory *elenent") due to the • diet’irbances '* (de- 
partures from •normal** conditions', 

2, First Differential Effects, 

The mathematical theory underlying the computa- 
tion of these linear approximations, or * diffarantial ef- 
fects * as they ore usually celled, was developed by G, A, 
Bliss while he was at the Aberdeen Proving Grcund during 
World War I, Some of his results were published In refer- 
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enc«s [l] md [s] of the bibliography at the end of this 
paper: others appear In references [s] » [d] , and [s]* 

Other investigators have provided modifications 
and improvements, both in the theory and in 
its ^plications. The most recent such contribution is 
a paper by S. J. McShane, entitled The Differentials of 
Certain Ftmotionals in Exterior Ballistics *, and soon to 
appear in the Duke Mathematical Journal. In this paper, 
which will hereafter be designated as [mcs], MCShane shows 
that, for most trajectories, the *norm* (or measure) of 
the disturbance as used by Bliss In. his proofs can be re- 
placed by a simpler, more natural "norm*. 

In order to make clear the distinction between 
the norm used by Bliss and that used by McShane, we shall 
quote below an excerpt from the introduction to the paper 
[mcs] { the manuscript of /Which Professor McShane very 
kindly loaned the writer in advance of its publication. 

To be specific, let us restrict our attention to 
the effect of wind on range. The 'normal* condition, un- 
der which trajectories are computed, is zero wind. The 
range is a functional of the wind. If the wind is sup- 
posed to be a continuous function of altitude alone, we 
are considering a real-valued functional on a space of con- 
tinuous functions. In this space we can introduce a 
familiar nom, the norm of a continuous function w(.) being 
defined to be the least upper bound of its absolute value. 
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Then the clasaical definition of a (first) differential can 
he jsrplied; it is a linear functional of v(,) which differs 
from the change of range due to w(,) by a quantity 

whose ratio to norm w(.) approaches zero with norm w{,). 

But the .existence of a differential in this sense has not 
■ eretcfore been established* To the best of my knowledge, 
tile < nly adequate mathematleel dis ..usslon of' differential 
ffects s due to G. A. Bliss* Bliss, however, does no? 

1^; ihe norm mentioned above. Instead, he defines the 
norm -n w(.) to be the, greater of the least upper bound of 
iwj and the leaitt upper bound of the absolute value of itr 
derivative with respect to altitude y; that is, he uses t'le 
norm of the Banach space C Instead Of the norm of the 
z,oace C, Now with certain disturbances this would be 
quite unobjectionable. For instance, the Coriolis force 
a project! e is a disturbance whose value and rates of 
change with respect to position and velocity are comput- 
a le, and the inclusion of the rates of change in the defi- 
nition of ~ho norm produces no annoyance. Bur. the wind is 
an experimentally determined function, and the experiment 
t .rnishes only the mean values between a finite number of 
successive altitudes. The experimental result may be ex- 
pe ted .0 differ from the exact wind by an error whose 
norm in the space C (i,a., maximum absojlute value) is 
sma.. but we have no assurance that the rate of cfange of 
wind with altitude does not at some places differ widely 
from the experimental estimates In order that the 
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mathematical treatment of differential affects shall be 
applicable to the experimental methods Involved, it is of 
great importance that the norm be determined h3r the magni- 
tude of |w| alone, and not depend on the derivative w’ . 

We ha'^e seen that disturbances can be classi- 
fied into two aggregates; for those of one class it is 
harmless to admit the derivatives in the definition cf the 
norm; the others, specifically wind, temperature, and den- 
sity, are to be regarded as functions of altitude alone, 
and their derivatives are to be excluded froFi the defini- 
tion of the norm. It is possible to avoid the annoyance 
of considering two separate desses of disturb'^nces, and at 
the same time t;» gain slightly in generality, by defining 
the norm of a disturbing |function to be the greatest of the 
maxima of the absolute values of the function itself ana of 

r 

its partial derivatives with respect j;o allf variables ex- 

f, ■» 

cept the alt%t’-de y, , . . „ , fhe norm being defined, 

with given 'normal’ oonditions and given dis- 
turbances, do departures from 'normal' condi- 

tions produce differentiable effects? We find that they 
do. . Is "the difference between differential ef- 

fect and actua.l change an infinitesimal of second order in 
the norm? V/e find that it is if the trajectory is every- 
where ascending or everywhere descending; otherwise' it may 
not be as small t>s second order," 
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3, Need For A New Theory. Svucmary of Results. 

For most purposes, tlie use of differential effects 
in place of actual effects of disturbances has produced 
quite satisfacto results* Recently, however, evidence 
has accumulated to support the belief that linear approxima- 
tions to effects of cer ain disturbances (e.g. , a wind; or 
the combined presence of wind and abnormal temperature) do 
not give sufficiently accurate estimates of the actual ef- 
fects. Because >4 the need of Improved estimates, and be- 
cause the literature contains no mathematical theory of 
second-or er effects of disturbances, Professor McShane 
suggested to the writer that he undertake the development 
of such a theory. 

'‘ The present paper consists of a portion of the re- 
sults of that investigation. It includes: (1) a proof of 
the existence of a second differential of the mapping de- 
fined by vhe dif .eren-^ '.al equations of motion of the mass- 
center of a proj ctll.:,; (2) a demonstration rhat this se ond 
differential is expressiblo in torms of the second variations 
of the coordinate', f the center of mass; (3) a derivation 
of the system of differential equations satisfied by the 
second variation^, a system which tas the same homogeneous 
terms as d- es the classical system of first variation equa- 
tions ; 4) a fcrmula for expressing a solution of ti)e dif- 

fareatial system or second 3ariat..ons in terms of double 
Stieltjes iitesrals, (5) a simple formula for calculating 



7 



the second verietiuns of traJecLorv eiemerite from the second 
varintions of the coordinates of the ioasfl-center ; and (o) a 
theori!!!. on second dif ferontlsls enalo<^us to ;..cJhane's main 
theore:\, (11.7 j of [l.co] , on first differentials. 

The nor." of disturbance employed In obtain in 3 the 
I'irrt J’ive results above is defined in terms of the absolute 
values of the disturbance functione and their partial de- 
rivatives of the first and second orders, v/lth respect to 
all the variables except the independent one, A somewhat 
sijapler norm Ir used in establishing? (d). 

Field s for ‘"urther Investigation. 

It has Ion' been known that (first) differential 
effects are expressible in terns of (single) Stleltjes In- 
tegrals, (By *'otieltjes Integral’*— both the single and the 
double varietie. — v.e mean the one often referred to as the 
"Aiemann-Stieltjes” integral, not the 'TLebesgue-stlelt Jes" 
integral.) This fact has made possible the use of a "bal- 
listic wlnd"--e. fictitious constant wind which has the same 
differential effect on range as does the actual range-wind 
Wx(y)» (A "ballistic density" and a **ballistic tempera- 
ture" are similarly defined.) The ballistic wind is com- 
puted ^ the field from observed wind-data by the use of 
so-called "wind weighting-factor functions" v/hlch have been 
pre-computed at the ballistics laboratory. iiach wind 
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Vifei>‘ 5 Ut lug- factor function is, by definition, the ratio of 
the differential effect, g(k), 0 '* k = 1, of a range-wind 
waich is zero at all levels above kY (Y being the altitude 
of the suiPj’Jlt of the trajectory) and +1 at all levels be~ 
low I'Y, to the differential effect g(l) of a range-wind 
which is one length-unit per time-unit at all levels* The 
graph of a weighting-factor function is known as a 
"weight lug-factor curve", and is plotted at the laboratory. 
Thus the fact that (first) differential effects ax‘e ex- 
pressible as (single) Stieltjes integrals has rendered 
possible a method whereby practically all the computations 
for differential effects can be performed by a group of 
people who can be situated far. from the combat theater, and 
the results of their calculations can be recorded in such a 
manner that applications may be made in the field quickly 
and by means of very simple computations. (?or a more de- 
tailed account of weighting factors and their' applications, 
see ]feferances [5], [6], and [7].) 

Now the results of the present investigation, 
which enable us to express second differential effects in 
terms of double stieltjes integrals, lead us naturally to 
conjecture that it is possible to devise a method of 
weighting-factor functions for second differential effects 
analogous to those for first differential .effects described 
above. But no attempr has been made to do this in the 
present paper. In fact, such an effort does not appear 
feasible until the formulas developed herein shall first 
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have been applied to the balliatlc aqusrions (differential 
equations of mc^^ion of the mass^oenter ) themselvea. This 
latter task is Itself one of some magnitude, as will become 
increasingly evident in later chapters. It has been post- 
poned as constituting an appropriate objective for a Inter 
investigation . 

5. Remarks on Notation. 

After experimenting with several notations* the 
writer concluded that, for the purpose of this paper, the 
most satisfactory system to use would be a slight modifica- 
tion of the notation employed in ^MoS] * Since the latter 
paper is not yet available in published form, it has seemed 
advisable (for the convenience of the reader) to quote sev- 
eral passages from the manuscript form of flJcs] • These 
will be clearly indicated by quotation m&rks, and grateful 
acknowledgment is hereby made for Professor llcShane’s per- 
mission to make these quotations. 

We shall use the summation convention with regard 
to repeated indices throughout the paper, except in certain 
parts of Chapter 71 where the departure from this notation 
is clearly indicated. 
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Chapter II 

THE EQUATIONS OF MOTION. A MORB OaSNERAL 
DlggERZHTlAL SYSTEM, 



6o The SQuatlons of Motion, 

In this section v;e shall develop the differential 
equations of motion of the mas s« center .of a projectile, sub- 
ject to certain simplifying assumptions, in the form in 
which they are usually encountered in ballistic literature. 
For the notational reasons indicated in ^ 5, we quote the 
following passage from [Mcs] • 

liet us choose a rectangular coordinate system in 
which the direction of the earth's gravity field at (0,0,0) 
is along the negative y-axls. The system is supposed fixed 
v/lth respect to the earth. In order to avbld I nes sential 
annoyances, we shell select some number B and decide to con- 
fine our attention to trajectories along which the in- 
equalities 

(6.1) lx|<B, ly|<B, U|<b; 1vxI<B, |vy|<B, |Vz|<B, 

are satisfied and whose time of transversal is less than B* 
The set of all points satisfying (6,1) will be named R, 

This restriction is not unrealistic; at any particular 
epoch there is an upper bound to attainable velocities, and 
if trajectories depart more than a few hundred miles from 
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the earth's surface problems of ballistics blend ^nto those 
of orbit computation. The components of velocity of the 
center of mass of the projectile (which for brevity we 
shall call the components of velocity of the projectile) 
will be denoted by (v^* Vy, v^) , and the length of the ve- 
locity vector by v. The components of wind- velocity will 
be denoted by (w^, Wy, Wz), It is convenient, though not 
absolutely essential, to assume that wy is identically zero. 
The velocity of the projectile with respect to the air then 
has components 

(6,2) Ux “ '^y “ ^y~''''y ~ ^y» '^z ” V2-W2 « 

The length of the vector (ux»Uy,U2) we denote by u, V/e 
a s sume that the components c f v;l nd -veloc ity , the air - 
density ^ , and the absolute temperature S ere f uno ,-cas 
of position only , independent of time . We shall assume 
that for a projectile of given shape, size, and mess, mov- 
ing through air of given chemical composition, the princi- 
pal aerodynamic force has direction opposite to the vector 
(uj,Uy,U2) euid magnitude muE, where m is the mass of the 
projectile and E is a function of u, ^ and anc is 
continuous together with its partial deri.etives of first 
and second order.* If the gravitational accelert ti^on were 
of constant magnitude and direction throughout it would 



* For a more detailed discussion of the nature of the func- 
tion E, see references [s] , [o] , and [?]. 
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produce a force (0, -niK» 0) on the projectile. Thus the 
force on the projectile has courjonents approximutely equal 
to “DxEuxf -niSuy-rag, ~n.dJ 2 . All the othei* forces, whether 
uv rodynattiic ■ r due to -ohe ch^iige in gravity-field or 
Coriolis forces, v/e lump tC( t'rer and designate oy the sym- 
bol (max, may, ^z) • I'or pre . t-i, . pux'poses it is immaterial 
how these forces can be calcul; ied; it is enough tlust they 
can be found, and that they ax*; ordinarily considerably 
smaller than ’drag’ mEu and gi*. rity mg. We can now wite 
the equations of motion in the form 





dt 


- Euy - 


g ay, 


ct 


“ Vy, 


(d,3) ^ 


dV;c 

1 ' 


- EUx -i- 


®XS 


c.x 

U u 


= Vx, 


1 

\ 


dvz ^ 
, ^t 


- 2U2 + 


aj.,, 


d2 

(ft 


= Vj,. 



V 



"’Henceforth we shell acsume that there is a oOil- 
• iv number ^ such that e voiyvn/faere in R the inequality 




holds. Then at every value of t such that Uy » Vy = 0, 




■g 







hove 
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This implies that whenever Vy vanishes y has a proper 
maximujQ, Then y can have no relative minimum. Between 
two proper maxima a minimum must occur; therefore there can- 
not be two times at which Vy vanishes. If on a trajec- 
tory there is a time ts» at which Vy » 0, the point 
(x(ts)* y(ts)» z(ts)) is called the summit. For t < ts 
vje have Vy > 0; the corresponding arc of trajectory is the 
ascending branch . For t >■ tg we have vy 0; the corres- 
ponding arc is the descending branch . Artillery trajecto- 
ries, for fire at ground targets, ordinarily have summits. 
Bomb trajectories in level or dive bombing have the descend- 
ing branch only; anti-aircraft trajectories usually consist 
of an ascending branch caily. 

"The assumption (6.4) is a very mild one for con- 
ventional trajectories. However, by adopting it we exclude 
projectiles . , . which are capable of sustained level 
flight." 

7 . T he Equations of Motion as a Special Case of a More 
General Differential System. 

The system (6,3) is only one of several ways in 
which the equations of motion of the mass-center can be 
written. Other useful forras are obtained by using alti- 
tude or slope as independent variable, instead of time. 

In each of these forms, hov/ever, the right members of the 
equations are rather complicated expressions involviiig the 



coiiipound junction £, But all these forms of the equation^ 
raotion hive in common the property of being a special 
'I'tce of a dif lerential system which can be written quite 
simply 1 ^ the notation introduced in [mcs] . Our investi- 
gation, t i-refore, will take the form of a study of this 
more general differential system, not only because of the 
greater econo i; of notation but also because the theory 
v;hich we .hal^ c.evelop will then be applicable to several 
forms of the equations of motion instead of being confined 
oaly to the system (6,3) » For this reason, we again quot : 
from [l';cSj , this time from | 4 of that paper. 

The variable v/hich is chosen as independent var- 
l. 'jle in the equations of motion shall be called s , and 
shall be assuiiiid to vary over some interval (sq,sg)o If h 
!r oime, this corresponds to assuming a bound on the timei 
vf flight of a ,1 projectile, considered; if, say y or x i ; 

' V. n 1 as indSi-endent ''a'-aole, it follows from the ass'^^ 
i' t I't the orajectori o lie in a bounded part of sp ce 
'ri.-. lepende.i 'ar^:. bl„s we shall q^,. . ,,q^ , and ch - 
i-'J.-.- (q^,. q*^) we shall abbreviate to q when Ul •: 

Sup-r.-.c j.pts aff X d directly to letters are 
:r;nr'.'?r 2 ve ; expone ts, v/hen we have occasion tr u, . 

-.’T ’ 11 be pf ixed to p ..renthe ses. Thus (q®)^ wool ? 

■* nqua.-f- of the thi 'd component of the n-uple q , 

■i'' - - ' r-1 e.xt V ria les ra. 7 (- over some open set Ro in n- 

- On the »t (si, g) Ro , consisting of all 
I r. . _.,njes (s,q) ’'■’ith sq ^ s = sg and q in Rq , certain 
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functions ^ (s,q), . . . , ® (spq) are defined; these 

will always be assumed continuous except perhaps at a finite 
number of values of s, at which they may have simple Jump 
discontinuities. We shall also assume that all the ^ 
considered are less in absolute value than some uniform 
bound B; this corresponds to deciding in advance that we 
shall only consider winds, velocities, etc,, which do not 
exceed a chosen bound. On the set in m-space defined by 
|r®^l “ B ( = l,...,m) there are n continuous functions 

F^(r F’^(r). We seek solutions of the equations 

(7.1) dq^/ds - F^( <f> (s,q(s))) 

m Fi( ipl(s, q^ , • . • , q*^ ),...,^p^(s, q^ , • • • , q^ ) ) , 

( i“l » . . • »n ) 



with initial values 

(7,2) q^(si) « qo (i»l,,..,n) 

lying in Ro and with q(s) (<l^(s) » •• • »<l^(s) ) 3 

lying in R© for sq * s - S2. 



*The equations (7,1) and the initial conditions 

(7.2) are determined when we specify the system of numbers 
and functions 

(7.3) (qj,...,qj, . . . ,f“) . 

"By well-known theorems [see, e.g., E. J. 
McShane: Integration , (Princeton Univ, Press), p, 344.] on 

differential equations, for each system (7,3) with the 
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h 

properties previous^ specified, equations (7.1), (7,2) have 
e solution defined on sorae sub-interval (si, 33 ) of 
(si, 83 ) and l 3 '’ing in Hq, But we cannot guarantee that 
the solution is unique. Syste'cs lacking unique solutions 
are of no interest to us, so from the aggregate of all sets 

(7.3) we select the subset (j) consisting of all 
(qo, ,F) vjhich furnish unique solutions.’* 

Since on the set (|) the solutions may not be 
continuable to s = sg, we select from (|) a subset (Jo de- 
fined as follov/s. 

(7.4) Definition : A "point" (qQ, , F) belongs to (Jo if 

it belongs to (J and the solution q(s) of (7.1), 
(7,2) is defined and in Rq for si * s = sg. 

That is, the equations (7.1), (7,2) map (Jo uniquely into 
the space On of n-uples of real-valued functions defined 
and continuous on the interval (si,sg), the "image" of 
(qo» f » F) being q( . ) = { q(s) | si = s = S 2 ^ , With 
reasonable noi-ms in (Jq and Cn, it can be shown that this 
mappin- is continuous on (Jq. In fact, this conclusion 
is established in Theorea (4.9) of.[r..cs] for the norms N 
and Nq defined below. 

("7 *5) Def inlt ion : Using the summation convention, v/e 

write Iql = (q^q^)^*^^. 



(7.6) Definition ; In the space Cq the norm will he the 
customary one, N(q(.)), given by 

N(q(.) ) * sup { lq(s) j | si - s = sg } 

(7.7) Definition ; In (fo the norm, which we shall 

designate No(qo» if i F) , will be the greatest 
of the three numbers ; Iq©!* 

sup { I If (s,q)l j Si » 3 » 32, q in Ro } *» and 

sup I |F(r)l I lr*^l < B, » l,...,m] . [By 

1 If I we mean ( if ^ similarly, 

lF| » j 

McShane's continuity theoreiUj on the mapping of ({iq into 
Cn» is as follows: 

(7.8) Theorem : When the spaces (J and Cn are metrized 

in accordance with the respective norms No and N, 
the set (j>o is an open subset of 0, and the 
mapping of (|)q into Cq defined by the equations 
(7,1), (7,2) is continuous. 

In the opening paragraph of this section, it was 
asserted that the system (6,3) is a special case of a more 
general differential system. That more general system 
consists, as the reader may already have inferred, of the 
equations (7,1) and (7,2), To see this, we identify the 
independent variable s v/ith t and the dependent variables 
qi with X, y, z, Vj, Vy, V 2 — -though not always in the order 
named. The functions we taHe to be: t, x, y, z, vx, 
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vy. ’?z» wx, W 2 , ^ , Q, ax, ay, a^. Thus, for the system 
(6,5), we have n ■ 6 and in ■ 14. Given a law of air re- 
sistance (i.e., a particular ’’drag function”), the function 
H can he calculated from the values of the 14 functions 

, and hence the right members of (6,3) are the func- 
tions with the r^ Replaced hy the 14 functions , 

which in turn are functions of the 6 dependent variables qi 
and the independent variable t. 

Theoretically there is a one-to-one correspondence 
between projectile-types and dreg functions. But the de- 
termination and tabulation of a drag function is an expen- 
sive process; not only are costly experiments (firing pro- 
grams) necessary, but extensive calculations are required 
to reduce the data obtained from the experiments. For- 
tunately, the function E has as a factor not only the drag 
function but also the reciprocal of an empirical quantity 
called the ballistic coefficient ; by altering the value of 
the ballistic coefficient when changing from one projectile- 
type to another, sufficiently similar to the first, ballls- 
ticians have found it possible to use the same drag function 
for more than one projectile-type. The functions F^ are 
always the same unless the drag function is changed; since 
it is rarely desired to calculate even first differential 
effects due to changes in drag function alone, we have not 
felt it worthwhile to complicate the notation in our theory 
of second differential effects by including terias which 
arise from changing the functions F^. 
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Accord lngly » the only ^^ disturbances ** which we 
shall consider in the remaining chapters are those which 
consist of changes In Initial conditions , or of changes In 
the functions ^ both . A disturbance which con- 

sists of changes both in the initial conditions end In the 
functions ip we shall denote by the symbol ( qo*^^ ) •' 

In accordance with this notation. ( A q.ot 0) will denote a 
disturbance which consists of changes in initial conditions 
only; ( 0, ) will mean a disturbance which consists only 

of changes in the functions . Therefore we shall 

suppress the symbols F and A F from such notations as 
(qo» 4 =^ » F) or No( A Qo. , a f). 
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Chapter III . 

THE (FIRST) VARIATION BQUATIOMa; CONTINUITY 
AND DIFFERENTIABILITY THEOREMS OF BLISS. 

6 L i ■ 

8. The (First) Variation gguatlbnsr 

At the beginning of $ 7 we made certain assump- 
tions regarding the continuity of the functions 

,m) and (1 ■ l,...,n)«. For our discus- 
sion of first and second differentials, of the mapping de- 

♦ 

fined by (7.1) and (7,2), we shall need additional assump- 
tions concerning the differentiability of the functions 
and pi. Those additional properties are described 
in the following 

(8.1) Definition ; Let (j^gg be that subset of (|)q which 
has the properties 

(8,li) the functions are twice continuously differ- 
entiable on the set of r with 
|r » B, •- l,...,m; and 

(8.1ii) for all s in (sijSg) and all q in R^, the functions 
tp ^ have bounded partial derivatives of the first 
and second ox-ders with respect to the qi, continuous 
except perhaps at a finite number of values of s. 

It should C; pointed out hei*e that the assumptions con- 
cerning seoG,. . derivatives contained in (8,li) and (8,lii) 
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are unnecessarily restrictive for the theory of first dif- 
ferential effects. In fact, McShane uses in [mcs] the 
subset ^3 of 4>o defined by altering definition (8,1) in 
such a way that the word "twice* is omitted from property 
(8,li) and the words "first and second orders" in property 
(8,liil are replaced by the words "first order". Since 
evidently (^33 c <Jg, the results in [mcS^ are valid on 
the set (Jgg. ] 

We wish to emphasize that in the definition (8.1), 
as in the analogous definition for (Js in [Mcs 3 » the var- 
iable s chosen as independent variable is distinguished 
from all the others, i.e., from the variables For in 

(8,lii) the partial derivatives of the functions with 

respect to £ are not mentioned. This is the reason for 
using the letter s in the subscripts on the letter (^. 

For example, the set (^tt (i,Q,» <tse for a - t) is not the 
same set as 4yy (the set (Jsa when y is chosen as inde- 
pendent variable ) . 

Now, let (qo» ^ ) and (5o» if ) points of 
4s SI and let 

C ^±^0 ^ \ 

(y.z) J 
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cf - if 



7or each real nunoer T ^ Q m f » 1, let 

s { • , t I ( s f *7 ) ) 

be the solution of the corref^onding differential system 



(8.5) 



d ejlj.*' 1^ ^ 

Is ~ 



s “f.- + 






i.e,, we assume that, for each t' , 0 = 1, 



(8.4) 



1 ds 






Cur 

0 



Let 



(9.S-J 



CS) 3 a? CS^o) 

Cs> - 




and, for br-evlty, let 
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+ 't A^tf Cs, Q.Cs^'tr))) 



M < 






[of r d^ . . .^ >T^ ; 

J = 1, . . . , >1 .] 



Z’* • 

if 



with siallar abbreviations for F ’'„ /a , ^ ^ k 

r«r<^ ’ 

A^ this notation, the result 

of differentiating each member of (4) with respect to 'f , 
and then interchanging the order of differentiation, may be 
v/ritten 















fl.7) 



•r 



Thus, the functions 

(8.8) s ri^Cs- ‘f. 



satisfy 



(8.9) 



( Jtfi'- , f'’ {°J 
) TT 



^ P 
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